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Abstract 

We investigate the prospects for indirect detection of fermion WIMPless dark matter at the 
neutrino telescopes IceCube and DeepCore. The dark matter annihilating in the Sun is a hidden 
sector Majorana fermion that couples through Yukawa couplings to a connector particle and a 
visible sector particle, and it exhibits only spin-dependent scattering with nuclei via couplings 
to first generation quarks. We consider cases where the annihilation products are taus, staus, 
or sneutrinos of the three generations. To evaluate the muon fluxes incident at the detector, we 
propagate the neutrino spectra through the solar medium and to the Earth and account for the 
effects of neutrino oscillations, energy losses due to neutral- and charged-current interactions, 
and tau regeneration. We find that for the stau and sneutrino channels, a 5 yr 3u detection of 
dark matter lighter than about 300 GeV is possible at IceCube for large Yukawa couplings or 
for dark matter and connector particles with similar masses. The tau channel offers far better 
detection prospects. However, due to its lower energy threshold and better muon background 
rejection capability, DeepCore is able to detect signals in all annihilation channels and for a 
wider range of dark matter masses. 



1 



1 Introduction 



There has long been interest in the use of neutrino experiments for indirect dark matter (DM) 
detection, with a particular focus on neutrino signals from DM annihilation in the Sun. The basic 
idea is that DM particles can be captured by the Sun if they lose sufficient kinetic energy through 
elastic scattering from solar nuclei. This leads to an enhanced DM number density in the solar 
core, where DM can annihilate to standard model (SM) products. These products in turn decay, 
producing a neutrino flux which can be detected by Earth-based detectors. Because the Sun is 
(for most models under consideration) in equilibrium, the DM capture rate is equal to twice the 
annihilation rate. Thus, the neutrino flux on Earth is determined by the capture rate, which in 
turn is determined (up to a few 0(1) factors related to solar physics) by the DM mass m x and 
the DM-nucleon scattering cross section a. This is one of the key advantages of neutrino-based 
indirect searches for DM: the event rate is directly related to the scattering cross section, without 
many of the astrophysical uncertainties which attend other indirect detection strategies. Moreover, 
it is in principle possible to directly compare search results at neutrino experiments with those at 
direct-detection experiments, which also measure the DM-nucleon scattering cross section; for a 
recent discussion, see pp. This allows one to corroborate a signal at one type of experiment with a 
signal at a very different type of experiment, which is desirable. 

The downside, however, is that for much of the (m x , (j)-parameter space, sensitivities for cur- 
rent and future neutrino detectors have already been surpassed by direct-detection bounds, im- 
plying that neutrino detectors may have difficulty providing new input to DM studies. However, 
there are two scenarios where neutrino experiments are expected to shine. One is at low-mass 
(m x ~ 4 — 10 GeV). In this mass range, bounds from direct-detection experiments tend to become 
significantly worse because the nuclear recoil energies often fall below the experimental threshold. 
Experiments such as Super-Kamiokande (Super-K) and other proposed water Cherenkov or liquid 
scintillator-type detectors can provide the best bounds in this range |2], and could be sensitive to 
spin-independent (SI) scattering cross sections asi ~ 10 -5 pb. DM in this range can potentially 
explain the annual modulation signal seen by the DAMA experiment ([3], though see also [3j), as 
well as unexplained events recently reported by the CDMS and CoGeNT Collaborations [5JE1E]- 
There is thus considerable interest to cross-check these results with neutrino detectors. 

Another scenario for which neutrino experiments can excel is for models where the DM-nucleon 
scattering is largely spin-dependent (SD), rather than SI [8]. The reason is that DM can then be 
captured by SD scattering from hydrogen in the Sun, leading to significant neutrino fluxes. Current 
direct-detection experiments provide much less sensitivity to SD scattering and the best bounds are 
set by indirect detection experiments, even at large m x . This has been the focus of large neutrino 
detectors, such as AMANDA or IceCube. For example, with 1800 live days of running, IceCube 
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could be sensitive to 0(100 GeV) DM with a S D ~ 10~ 5 pb [9]. Although the relatively high 
threshold energy of these detectors limits their utility in studying the low-mass DM, their large size 
makes them ideal for studying higher-mass DM (m x > 100 GeV) with ctsd 3> °~si- 

Direct detection is much more sensitive to o~si due to the possibility of coherent scattering 
in the nucleus. Direct-detection experiments, including CDMS and XENON100, currently have a 
sensitivity to osi in the O(100 GeV) mass range of about 3 orders of magnitude [10] greater 
than the IceCube projected 1800 d sensitivity to o~sd- But for models with a si <C ctsdi comple- 
mentary coverage from neutrino detectors becomes important. This complementary coverage is 
especially important if one goes beyond models of weakly interacting massive particles (WIMPs). 
The recently proposed WIMPless model of dark matter [T2] is notable because it provides a 
very generic hidden sector DM candidate which naturally has about the right relic density to match 
astronomical observations, regardless of the mass of the DM particle or the details of the hidden 
sector. This versatility suggests the possibility of specific WIMPless models with almost exclusively 
SD scattering, for which the best detection prospects would lie at neutrino detectors. 

We present a model of WIMPless dark matter where the DM candidate is a Majorana fermion 
in the hidden sector. At tree level, it will exhibit only SD scattering with SM nuclei. We find that 
significant neutrino fluxes on Earth arise from models where DM annihilates either to tau, stau, 
or sneutrino pairs. We consider the detection prospects for this type of model at DeepCore and 
IceCube. 

The paper is organized as follows: In Sec. 2, we review WIMPless dark matter and describe the 
properties of the fermion candidate DM particle. In Sec. 3, we obtain the neutrino spectra arising 
from the decay of various annihilation products. In Sec. 4, we discuss the predicted event rates at 
DeepCore and IceCube, and we present our conclusions in Sec. 5. 

2 WIMPless Dark Matter 

Here we briefly review the WIMPless DM scenario |11| . This scenario is closely related to gauge- 
mediated supersymmetry breaking (GMSB) |13| . Our theory consists of a supersymmetry (SUSY)- 
breaking sector, the visible minimal supersymmetric standard model (MSSM) sector and a hidden 
sector, plus a connector sector which we introduce in the next subsection. We assume there is 
at least one chiral superfield S in the SUSY-breaking sector, and there are messenger fields <£, 
$ between the SUSY-breaking and MSSM sectors and <J> X , $ x between the SUSY-breaking and 
hidden sector. S is coupled to the messenger fields through Yukawa couplings in the superpotential 
W = \®S<Z> + \ X $ X S$ X . The effect of GMSB on the MSSM sector is well-known: the chiral 
field S acquires a vacuum expectation value (S) = M + 9 2 F. The Yukawa couplings generate 
messenger mass terms of order m mess ~ AM and messenger mass splittings of order F mess = XF. 
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Once the heavy gauge messengers are integrated out, the low-energy effective theory has a new 
SUSY-breaking soft scale m so ft which is generated by diagrams with the messengers running in the 
loops, 

9 2 ( F mess \ _ g 2 F 
m *oft~ {A%)2 [ mmes J ~ (47r) 2 M > W 

where g is the largest relevant gauge coupling. The effect in the hidden sector is qualitatively the 
same, and we assume that some unbroken symmetry stabilizes a particle in the hidden sector at 
the soft SUSY-breaking scale, whatever it may be. The soft SUSY-breaking scale in the hidden 
sector is then given by 



2/tt, \ „2 

(2) 
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x (4tt) 2 \m messx J (4vr) 2 M' 
where F messx = X X F and m messx ~ A X M. Because the ratio F/M is determined by the dynamics 
of the SUSY-breaking sector, it appears in the same way in the soft scale in both the hidden sector 
and MSSM sector. We thus find 
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The ratio g 4 /m 2 sets the annihilation cross section for a particle of mass of order m through 
gauge interactions of strength g, and the cross section in turn determines the thermal relic density, 
(fi oc (<Tf) _1 [14J). To obtain a thermal relic density which matches astronomical observations, 
a DM candidate would need an annihilation cross section (av) ~ 1 pb. The "WIMP miracle" 
is the amazing fact that, when calculated at the electroweak scale, the ratio gw ea k/ m weak yields 
an annihilation cross section which is in fact close to 1 pb. But we see from Eq. (|3|) that the 
annihilation cross section for a stable particle at the hidden sector soft SUSY-breaking scale will be 
approximately the same; this DM candidate will thus automatically have approximately the right 
relic density. Most importantly, the scale m x of the DM is a free parameter, thus motivating DM 
searches over a wide range of masses. 

2.1 Fermionic Model 

In the WIMPless scenario, DM is stabilized by a hidden sector symmetry (which may be discrete). 
We let Xl^r be chiral supermultiplets in the hidden sector, with a fermionic mass eigenstate of 
Xl,r (henceforth denoted as X) being the lightest particle charged under the stabilizing symmetry. 
WIMPless DM interacts with the SM through Yukawa couplings, which permit the exchange of 
exotic "connector" particles. Consider the interaction superpotential, 



W = \ L iX L Y L Q Li + \ Ri X R Y R Q Ri + m Y Y L Y R 
\' Lj X L Y^-L Lj + 

A Rj 



+\' Lj X L Y l L ep -L Lj + X'^XnY^Lzj + m Y u P .Y l L ^Y l R ^ , (4) 
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where Qu,m are MSSM chiral quark multiplets, Lhi,Ri are MSSM chiral lepton multiplets, and 
Yl,r and Y l j^ R ' are the exotic connector particles which are charged under both the MSSM and the 
hidden sector symmetry. Gauge-invariance requires the Yj, r to be chiral under the MSSM, so their 
fermion components behave like exotic fourth generation quarks (similarly, the fermions of Y^ R ' are 
exotic fourth generation leptons). The mass terms my, m Y iep. are thus set by electroweak symmetry 
breaking, and are determined by the Higgs vacuum expectation value and Yukawa couplings. 

The mass of the exotic quarks is constrained by perturbativity and precision electroweak data 
to the range my < 600 GeV [15j HSJ H7] . Since the DM particle X is the lightest particle charged 
under the stabilizing hidden sector symmetry, it is constrained to be lighter than any of the exotic 
connectors. 

Since the exotic squarks Yl r can get mass terms which are decoupled from electroweak sym- 
metry breaking, there is no relevant upper bound on . But can be bounded from below 
by direct searches at colliders. The exotic squarks can be pair produced through QCD processes, 
with each decaying via Yl^r — > X + u, d. The signal at a hadron collider would be exclusive dijet 
production with missing transverse energy, which is the same signature as that of leptoquark pair 
production, with decay to a quark and neutrino. A CDF Collaboration search for this signature [18] 
places a bound on the leptoquark mass at m^g > 187 GeV. Since, in our case, the missing energy 
arises from a relatively heavy-DM particle X (as opposed to a nearly massless neutrino), the lower 
bound on the mass of the exotic squark will be even weaker, and for our analysis it will not be 
constraining. 

The DM candidate X can be either a Dirac or Majorana fermion. We thus see that DM-nucleon 
scattering can arise from s- or u-channel exchange of the Y multiplets, as represented in Fig. 1. 
After a Fierz transformation, we can write the DM-parton scattering matrix element in terms of 

H 

Figure 1: The u- and s-channel Feynman diagrams for Xq scattering. 
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an effective operator, 
O 



\ 2 



2{m] 



We assume from here on that the scalars from the Yl and Yr multiplets have degenerate mass; 
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rriy^ = rriy- If we only retain terms which provide a velocity-independent contribution to 



DM-nucleon scattering, we can write this operator as 
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We find that if X is a Dirac fermion, then 



<?SD 



4 QA.i-{m\ 



m 



x) 2 



Y J {\l l + \ 2 Ri ){ZB p l +{A-Z)Bl) 



1 



m,. 



J+l 



4 47r(m| 



2 \2 



J 



(7) 



where the 1/4 factor accounts for the fact that the DM is not its own antiparticle, so only the s- 
or n-channel will contribute to Xq — > Xq scattering. m r = rnxMpj / (mx + Mjv) is the reduced 
mass of the DM- nucleus system. The vector- vector part of the effective operator generates 057, 
while the axial-axial part generates &sd- The B^f' n ^ parametrize the quark content of the nucleon, 
and the A^'™^ its quark spin content. The (S p ^) are the expectation values of the spin content of 
the proton (neutron) group in the nucleus. The numerical values of corresponding to structure 
functions including next-to-next-to- leading order QCD corrections are |19j . 



A? = 0.78 ± 0.02 . 



-0.48 ±0.02. 



A? 



-0.15 ±0.02. 



(8) 



As noted in Ref. [20J, uncertainties in quark spin contents can lead to substantial variations in as d- 
If the DM is Dirac, then a asD which is potentially detectable at DeepCore would already 
be ruled out by direct-detection constraints on a si- We therefore focus on the case where X is a 
Majorana fermion which is stabilized by a Z2 symmetry (so CXC = X). The vector- vector part 
of the effective operator is then identically zero and the scattering cross section is necessarily SD 
(a si = 0), and we find 

2 
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since now both s and u channels contribute. Note that we have only included the velocity- 
independent terms, and only at tree level. Since the bounds on a si are ~ 4 orders of magnitude 
stronger than the expected IceCube bounds on c~sDi one might wonder if this model could be probed 
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by direct-detection experiments via velocity-dependent or loop-induced spin-independent effective 
operators. One would expect velocity-dependent operators to be suppressed by at least 10 -6 . Also, 
loop-induced SI effective operators would be suppressed by (a/47r) 2 ~ 10~ 6 . Thus, with further 
improvements, direct-detection probes of a si induced by one- loop diagrams may be comparable 
to probes of ctsd at tree level by IceCube. This is an interesting subject for future study, but is 
beyond the scope of the present work. 

In our model, DM can have axial- vector couplings to the light and/or the heavy quarks. Scalar 
DM coupling to heavy quarks was studied in detail in [21], but no such study has been performed 
in the case of axial-vector coupling. Instead, we consider the case where DM couples to light 
quarks; if DM couples to more than one quark generation, then there is a potential for dangerous 
flavor-changing neutral currents. As a simple solution to this constraint, we assume that the 
only DM-quark coupling is to the first generation. We take A = Xl u = ^Ru = ^Ld = ^Rd, and 
Als = ^Rs = 0. 

2.2 Annihilation 

In WIMPless DM models, the annihilation cross section to hidden sector particles is naturally 
(cr ann v) ~ 1 pb at the time of decoupling. From the form of the Yukawa couplings, we see that DM 
can also annihilate to MSSM matter multiplets. The annihilation cross section to MSSM particles 
should not be very much larger than 1 pb, or else annihilation will dilute the relic density and ruin 
the WIMPless miracle. But in order to get a significant flux at neutrino detectors, it is necessary 
for the annihilation branching fraction to MSSM particles to be comparable to unity. 

Since we are considering the case where the DM is Major ana, annihilation to light fermions is 
chirality/p-wave suppressed. Also, in order to obtain SD scattering without inducing dangerous 
flavor-changing neutral currents, we assume that the only quark multiplets which the DM couples 
to are the first generation multiplets. So DM annihilation to quarks will be suppressed; if the DM 
annihilation to the hidden sector is not suppressed at current times, then the branching fraction 
for DM in the Sun to annihilate to quarks will be very small. One caveat to this statement is that 
if the light hidden sector particles are also fermions, then annihilation to the hidden sector is also 
chirality/p-wave suppressed. This does not significantly suppress the annihilation cross section at 
decoupling, but may suppress the annihilation cross section to hidden sector particles in the Sun 
by enough to allow the branching fraction to quarks to be sizable. In this case, DM annihilation 
to quarks in the Sun can be a significant source of neutrinos observed at IceCube. But if there are 
light vectors or scalars (perhaps Goldstone bosons) in the hidden sector, then annihilation to these 
products is not suppressed, and the cross section for hidden sector annihilation would be expected 
to be ~ 1 pb. To get an annihilation cross section to SM particles of comparable magnitude, one 
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Figure 2: Feynman diagrams for DM annihilation into r r + and into fifa. 



must either consider squark final states, or make use of Yukawa couplings to SM leptons. Then, 
the only relevant MSSM annihilation products are r, II r and the first generation squarks ul,r and 
d-L R- These annihilation cross sections are given by 
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and the corresponding diagrams are shown in Fig. 2. We set \' T = \' L 



X' Rt and A- 
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A' 



Ri • 



DM annihilation to squarks requires the squark masses to be lighter than the DM mass, which 
in turn is lighter than the connector, which is bounded by precision electroweak data and pertur- 
bativity to be lighter than 600 GeV. Generally though, in GMSB-inspired models the squarks are 
heavier than the sleptons due to the hierarchy of the gauge couplings. Thus, such light squark 
masses are not natural. Besides, the cascade products of heavy squarks and the resulting neutrino 
spectra are extremely model dependent, being determined by the features of the sparticle spec- 
trum. For these reasons we decide to avoid the treatment of annihilation to squarks and focus on 
the channels dominated by leptonic Yukawa couplings. 



3 Annihilation Products and Propagation 
3.1 DM Annihilation in the Sun 

The time dependence of the number of DM particles in the Sun is determined by the balance of 
capture and annihilation. A massive particle can be gravitationally captured by the Sun from the 
galactic halo at a rate Cq. On the other hand, annihilation at the center of the Sun will decrease 
the number of particles at a rate Ta- Thus, the time evolution of the number of DM particles in 
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the Sun ./V is given by the solution to the differential equation [22J, 

d -^ = C & -C A N\ (11) 

where C A = 2F A /N 2 . 

Assuming that the number of collisions DM particles undergo inside the Sun during the Sun's 
lifetime is large enough for them to thermalize |23j , the number density of DM particles at a distance 
r from the solar core can be expressed in terms of the DM mass mx, the temperature of the Sun 
T and the gravitational potential (f)(r) = 2npr 2 /(3M|,j): 



n(r) = n e"^ (r) ^ = n e 3M « 1 , (12) 

where no is the DM number density at the center of the Sun, p is the average density of the Sun 
and Mpi = 1.22 x 10 19 GeV is the Planck mass. C A is then a constant that depends only on p, T 
and the annihilation cross section (ctav) averaged over the velocity distribution in the limit v — > 0: 

fn(rf(a A v)d 3 r f SM^T \^ 

CA = (1 f \J3 \2 = \ aAV > • ( 13 ) 

(J n(r)arry \ mxp J 







N(t) = W^tanhf- ) , (1 i) 



Equation (jlip admits the solution, 



C A 

where r = 1/^CaCq is the characteristic time necessary to reach equilibrium. For t ~ 4.5 Gyr and 
{&ann.v) ~ 1 pb (as is the case for WIMPless dark matter) and ctsd^ 10 -6 pb (the range of interest 
for IceCube/DeepCore), it is known [22j[23] that t/r > 1 for the range 150 GeV < m x < 600 GeV 
considered here. Since T A = CaN 2 /2 = C Q tanh 2 (t/r)/2, one can easily see that when t » t, 
Ta — > T eg = Cq/2. 

This limit simplifies the calculations, as it relates the event rate to the DM-nucleon scattering 
cross section, thus bypassing all the astrophysical uncertainties related to the solar model. How- 
ever, the condition of equilibrium is not guaranteed. For example, in Ref. [23] it was recently 
shown that there are regions of mSUGRA parameter space for which the annihilation rate is far 
below the capture rate. Moreover, since DM can annihilate to hidden sector particles as well as 
MSSM particles, the annihilation rate relevant for neutrino detection is scaled by the branching 
fraction to MSSM decay products, Y^ SSM = T aBj? . In the figures presented in Sec. 4, we 
account for these uncertainties by scaling the muon event rates at the detector by the parameter 
£ = T¥ SSM /Y eq = Bp ISSM tanh 2 (t/r) . However, our discussion assumes £ = 1, as this choice al- 
lows us to draw quantitative conclusions. Event rates at conditions away from equilibrium can be 
obtained by simply choosing a lower value for £. 
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3.2 Capture Rate 

If WIMPless DM is Majorana, it has only SD scattering. This property presents the advantage of 
not being tied to the direct-detection bounds on the SI cross section. The corresponding capture 
rate is [221125]. 



where we have taken the local DM density to be 0.3 GeV/cm 3 , and the root-mean-square of the 
velocity dispersion in the halo to be 270 km/s. For SD scattering in the Sun, the only relevant 
nucleus with spin is hydrogen. 

3.3 Neutrino Spectra 

We limit our consideration to the cases of DM particles annihilating to taus, staus and sneutrinos of 
all three families. For the cases where DM annihilates to sparticles, we assume that the annihilation 
product is the next-to-next-to-lightest SUSY particle (NNLSP), and the lightest neutralino is the 
next-to-lightest SUSY particle (NLSP); in GMSB the lightest SUSY particle (LSP) is the gravitino. 
Such a choice is consistent with generic sparticle spectra originating from GMSB. 

In Figs. 3a, 3b and 3c, we show the neutrino energy spectra dN/dE v for an annihilation of 
WIMPless particles of mass mx = 150 GeV into taus, XX — > t~t + . The corresponding antineu- 
trino spectra are shown in Figs. 3d, 3e and 3f. Figures 3a and 3d show the spectra at the origin; 
Figs. 3b and 3e at the surface of the Sun and Figs. 3c and 3f at the detector. Similar spectra are 
obtained for annihilation into taus for the other two benchmark DM masses we consider in this 
paper: rax = 300 GeV and rax = 400 GeV. In Fig. 4 we show the neutrino and antineutrino 
spectra from annihilation into staus (XX — > t{"7j ) of a mx = 300 GeV DM particle. In Fig. 5, 
the neutrino and antineutrino spectra from mx = 400 GeV DM annihilation into VeLVeL, ^L^/iL 
and v t lVtL with equal branching fractions are shown. The spectral shapes for mx = 150, 400 GeV 
DM annihilation into staus are similar to those in Fig. 4, while those for mx = 150, 300 GeV 
DM annihilation into sneutrinos are similar to Fig. 5. 

We obtain the neutrino spectra by cascading the annihilation products with PYTHIA 6.4 [26] . 
In the cases of Figs. 4 and 5 we consider a typical GMSB spectrum, like the one given in Ref. [27] . We 
assume = 137 GeV and m^o = 94.5 GeV, corresponding to tan/3 ~ 10. The cascade proceeds 
primarily through fj — > Zfr~( + \ hence the higher number of u T s and v T s at the origin in 
Figs. 4a and 4d. 

In the case where the DM candidate annihilates primarily into sneutrinos (Fig. 5), we assumed 
for simplicity that the three sneutrino masses are degenerate (m^ = 111.5 GeV) and that the 
couplings to the three flavors are identical, so that the production spectra for the three neutrino 
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Figure 3: Neutrino spectra and fluxes for all three flavors per XX — > t~t + annihilation with 
mx = 150 GeV, a) at production, b) at the surface of the Sun, c) flux § v = N/AttR^ at Earth, 
d), e), f) v spectra and fluxes. 




Figure 4: Similar to Fig. 3 for XX -> f x with m x = 300 GeV. 

families are the same, v e L(jiL,TL) ~* ^Ve^r)- As m the previous case we have m^o = 94.5 GeV, 
corresponding to tan f3 ~ 10. Note that since WIMPless DM particles only annihilate leptonically 
into SM particles or the NNLSP, the neutrino spectra have the same shape whether X is heavier or 
lighter than the top quark. Notice that the spectra from tau decays are broader than in the stau 
and sneutrino cases. This is expected because of the influence of the neutralino mass on kinematics. 
As a matter of fact, the steep drop in the neutrino spectra arising from annihilation to sparticles 
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Figure 5: Similar to Figs. 3, 4 for XX 
m x = 400 GeV. 



r)^e(/i,r) with equal branching fractions and 



will have important consequences for detection at IceCube and DeepCore, as we discuss in Sec. 4. 

The propagation of neutrinos produced at the center of the Sun through the solar medium 
and to the Earth is detailed in the appendices. We take into account neutrino oscillations and 
the effects of neutral-current (NC) and charged-current (CC) interactions and tau regeneration. 
In Figs. 3b, 3e and corresponding panels of Figs. 4 and 5 we show the neutrino and antineutrino 
fluxes after propagation to the surface of the Sun. They present expected features; for example, an 
accumulation of events at lower energies is visible, due to energy losses of the neutrinos undergoing 
NC and CC scattering inside the Sun. Tau regeneration, too, has the effect of increasing the low- 
energy neutrino flux. Furthermore, we see that flavor mixing affects the and v T spectra but 
leaves the v e spectra unaltered, since #23 = 45°, #13 = 0° (Eq. [18] of Appendix A). Figures 3c, 
3f and corresponding panels of Figs. 4 and 5 show the fluxes at the Earth after our averaging 
procedure, described in Appendix A. We define the flux <]?„ = N ' j \Ai:R 2 SE ), with R$e ~ 1 AU being 
the Sun-Earth distance. Our results are in very good agreement with those of Ref. [28]. Notice 
that almost all modulation is washed out. 

The flux of high energy neutrinos produced by DM annihilation can be detected by neutrino 
telescopes like IceCube [29]. The IceCube detector at the South Pole is an array of 80 strings 
uniformly spaced 125 m from one another, each with 60 digital optical modules. It is deployed at a 
depth between 1450 m and 2450 m, and the instrumented volume of ice covers 1 km 3 . When muon 
neutrinos undergo CC scattering on nucleons in ice they produce a muon that propagates through 
the ice. These muons leave a trail of Cherenkov radiation that can be detected by the optical 
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modules. We are therefore interested in determining the muon and antimuon flux in and near the 
detector, taking into account possible energy losses due to the propagation in ice. DeepCore [30J, 
an extension of IceCube, is a denser (72 m spacing) array of six strings surrounding one of the 
central strings of IceCube, most of them placed below the "dust layer" at a depth of 2100 m. 
As detailed in Appendix B, DeepCore's denser strings lower the neutrino energy threshold with 
respect to IceCube and thus allow detection of neutrinos from annihilation of DM particles of a 
lower mass. Besides, DeepCore is designed to use the outer instrumented volume of IceCube as a 
veto on atmospheric muons from above the horizon to a level of one part in 10 6 , thus drastically 
reducing the background. Such properties make it a perfect instrument to test the WIMPless 
scenario in the case of Majorana DM. 

DeepCore's muon effective volume becomes insignificant below muon energies of 10 GeV |31j . 
At energies between 10 and 35 GeV DeepCore cannot provide directional information, since only 
one or two optical modules will be triggered. Very good directional information is needed to track 
the Sun throughout the year, so we fix DeepCore's energy threshold at E m i n = 35 GeV [32J. For 
IceCube we consider the very optimistic threshold E t h r = 100 GeV. 



4 Results 

In this section we investigate the prospects for DM detection at IceCube and DeepCore. We consider 
a 3cr-detection as 

"3cj-detection" o =g = 3 , (16) 
V N BG 

where is the observed number of muon events and Nbg is the number of muon events due to 
the atmospheric neutrino background |33| . which is discussed further in Appendix C. 

For the IceCube detector we divide the muon events into upward and contained events, fol- 
lowing the treatment of [34J. As we discuss in Appendix B, the former are due to upward going 
neutrinos interacting outside the detector volume, the latter to neutrinos that interact within the 
instrumented volume. To evaluate the event rate we track the Sun in the sky at different times of 
the year, and consider only events detected when the Sun is below the horizon. The distinction be- 
tween upward and contained events does not apply to DeepCore, which only sees contained events. 
For DeepCore, IceCube will veto atmospheric muons from above the horizon to a level of one part 
in 10 6 . Therefore, at DeepCore we consider the signal detected throughout the year. We consider 
a period of 5 yrs of observation. As anticipated in Sec. 3, we rescale the event rates obtained at 
equilibrium by £ = T^ SSM /T eq . 
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a ) tt, m x =150 GeV b ) rr, m x =300 GeV c ) tt, m x =400 GeV 




rrif (GeV) m i (GeV) m ? (GeV) 



Figure 6: Number of upward and contained events in 5 yrs of observation at IceCube for three 
values of mx as a function of m Y for the t~t + channel for several choices of A = \{i,,m{u,&)- 
Dotted black: A = 0.05; dashed blue: A = 0.1; solid thin red: A = 0.3; solid thick red: A = 0.5. 
The blue horizontal lines represent the number of events needed for a 3a discovery in 5 yrs of 
observation. The dashed brown horizontal lines represent gsd m pb. The yellow band represents 
values of A that yield the correct relic abundance. 

4.1 Event Rates: Tau channel 

In Fig. 6 we plot the upward and contained event rates at IceCube for our three benchmark DM 
masses for the XX — > t~t + channel, as a function of the scalar connector mass m-p and the Yukawa 
couplings. We assume that the masses of the up-type and down-type squarks Y are degenerate, 
and that A = \(L,R)(u,d) ■ I 11 Fig. 7 we plot the contained event rates at DeepCore for the same 
channel. The event rates [see Eqs. ([35]) . ([37]) and ([39]) of Appendix B] are evaluated for values of 
the Yukawa couplings which allow calculations in the perturbative regime. The values in pb of SD 
cross sections that correspond to the event rates are shown as horizontal dashed brown lines, and 
labeled on the plots. Notice the following features: 

• As we show in Appendix C, the atmospheric background to upward events in IceCube is 
~ 6.1 yr _1 , while the background to contained events is N^q 1 ~ 15.6 yr _1 . The background 
to contained events at DeepCore is NgQ — 2.5 yr _1 . The blue horizontal lines in Figs. 6 and 7 
indicate the number of signal events needed for a 3a detection in 5 yrs. We see that a 3cr-detection 
in upward events at IceCube is possible for reasonable values of A and niy- The contained event 
sample offers even better prospects. 



14 



a ) tt, m x =\50 GeV b ) rr, m x =300 GeV c ) rr, m x =400 GeV 




m f (GeV) m ? (GeV) m f (GeV) 

Figure 7: Similar to Fig. 6 for contained events at DeepCore. 

• In Fig. 7a we see that DeepCore significantly improves the prospects for observation of events 
for mx = 150 GeV. A 3<r detection in this channel can be obtained with Yukawa couplings of about 
0.3 without a need for resonant enhancement. If the Yukawa couplings are allowed to assume larger 
values, the prospects for signal detection become even more robust. 

• For DM masses well above E t \ ir the advantages of DeepCore in the r channel are less evident. 
As a matter of fact a 3a detection is more likely in the contained events at IceCube. For such high 
masses the propagator suppression of the SD cross section reduces the flux at the detector and 
makes the larger effective volume of IceCube advantageous. 

One may ask whether it is possible to relate the size of the Yukawa coupling A responsible for 
SD scattering to the size of the Yukawa couplings A T that appear in Eq. (jlpp . and are responsible 
for the annihilation to SM particles. This is possible under certain assumptions: 

• The yellow bands in Figs. 6 and 7 indicate the region with 0.1 pb < a xx ^ Tf v < 1 pb, 
assuming A' T = A and my iep . = m Y - For this range of cr xx ^ Tf v (assuming there are light bosons in 
the hidden sector), annihilation to MSSM particles has a large enough branching fraction to permit 
observation at neutrino telescopes, while not being so large as to dilute the relic abundance and 
thwart the WIMPless miracle. Note however, that this constraint is only relevant if DM annihilation 
to the hidden sector is not chirality/p-wave suppressed. 

• If X' T /X > 1, the yellow band shifts down to lower values. Such a situation may be desirable 
in the r~r + channel so that the iso-A curves pass through the yellow bands, i.e., the relic abun- 
dance and observable signals at neutrino telescopes can be simultaneously obtained for more of the 
parameter space. 

• As is clear from inspection of Eq. @ and the first of Eqs. (|10p . the capture rate depends on 
A and m Y while the annihilation cross section to taus depends on A^ and m Y ie P . • £ depends on the 
details of the solar model, particularly on the ratio T/p in Eq. (|13p . and on the ratios X' T /X and 
m Y /m Y i ep ,. A detailed analysis of these aspects is beyond the scope of this paper. 
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a ) f j f j , m z =300 GeV b ) f x f > , m x =400 GeV 




m 9 (GeV) m f (GeV) 



Figure 8: Similar to Fig. 6 for the Tj f± channel. IceCube is insensitive to this channel for 
nix = 150 GeV. The thick pink horizontal line represents the Super-K 90% C. L. upper bound 
on a sd- 



sl ) f , f , , m x =150 GeV b ) f 1 f 1 , m x =300 GeV c ) f l f 1 , m x =400 GeV 




nty (GeV) m ? (GeV) m f (GeV) 



Figure 9: Similar to Fig. 7 for the fi channel. The thick pink horizontal line represents the 
Super-K 90% C. L. upper bound on ctsd- 

4.2 Event Rates: Stau and sneutrino channels 

In Figs. 8 and 9 we compare the performances of IceCube and DeepCore, respectively, for the 
stau channel. Figs. 10 and 11 show our results for the sneutrino channel. The 90% C. L. upper 
bound placed by Super-Kamiokande [35] on asD is indicated on the plots by a thick horizontal 
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Figure 10: Similar to Fig. 6 for the vv channel. IceCube is insensitive to this channel for 
mx = 150 GeV, 300 GeV. The thick pink horizontal line represents the Super-K 90% C. L. up- 
per bound on (Jsd- In a), the solid blue horizontal line is coincident with the dashed brown line 
representing &sd = 10~ 2 pb. 
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Figure 11: Similar to Fig. 7 for the vv channel. The thick pink horizontal line represents the 
Super-K 90% C. L. upper bound on ctsd- 

pink line. For m x = 150 GeV this bound is afg x ~ 3.5 x 10~ 3 pb and for m x = 300 - 400 GeV, 
a Max _ g x 1Q -3 pb _ FoJ . the 

r t + channel these bounds are off scale and therefore not visible. 
The advantages of DeepCore in the observation of neutrino fluxes from low-mass DM are evident: 

• For our choice of spectra and energy thresholds, DM particles with a 150 GeV (300 GeV) mass 
do not produce observable events at IceCube in the two sparticle channels (sneutrino channel). At 
production, neutrino spectra for annihilation of mx = 150 GeV DM particles become negligible 
above 60 GeV for the stau channel and above 40 GeV for the sneutrino channel. Neutrino spectra 
originating from mx = 300 GeV DM annihilation into sneutrinos are negligible above 80 GeV. 
Thus, with a threshold energy E^. = 100 GeV no signal is expected for these channels at IceCube. 

• As Fig. 8a shows, even if observation of upward events at IceCube is possible in the stau channel 
for mx = 300 GeV, the region of parameter space allowed is very narrow. It is represented by the 
strip between the thin blue and thick pink horizontal lines. Similarly, Fig. 10a shows that a 
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3a discovery in upward events for mx = 400 GeV in the sneutrino channel is incompatible with 
Super-K data, as the thin blue line is above the thick pink line. Figs. 8c and 10b show that, as usual, 
the contained events offer better detection prospects at IceCube: assuming A = 0.5, 3a detection 
is obtained with a connector mass within 31% of mx = 300 GeV in the stau channel, and within 
24% of mx = 400 GeV in the sneutrino channel. 

The performance of DeepCore is much less sensitive to the annihilation channel than IceCube, 
because the lower energy threshold allows integration of much of the neutrino spectra independent 
of its shape or features. 

• Even for values of mx that do not allow observation of events at Icecube, DeepCore is able to 
observe events at 3a in 5 yrs for reasonable values of A and my. However, the region of parameter 
space for which DeepCore is sensitive to the mx = 150 GeV sparticle channels is very narrow, as 
can be seen in Figs. 9a and 11a. 

• The advantages of DeepCore over IceCube are not substantial for heavy-DM masses in the 
stau channel, and the detector fares just a little better in the heavy-DM sneutrino channel. Fig. 9c 
(11c) shows that 3a detection of mx = 400 GeV DM in the stau (sneutrino) channel can be 
obtained with a range of connector masses extended by a mere 3% (18%) with respect to IceCube. 

• The yellow bands represent the region of (my, A)-parameter space such that 0.1 pb < 
a ^, z ?v <1 pb, assuming that A' = A and m Y ie P . = my, 

It has been recently pointed out [36j that DeepCore should be able to identify contained "cas- 
cade" events, which originate from CC interactions of electron and tau neutrinos inside the detector, 
and from NC interactions of neutrinos of all three flavors. For such events the signal is enhanced 
with respect to the background at energies > 40 GeV, since the signal is predominantly due to the 
CC interactions. Specifically, at such energies the flux of atmospheric is from one to a few orders 
of magnitude larger than the flux of v e and v T [37]; thus, give the dominant contribution to the 
background through NC interactions, which are weaker than CC interactions. This could provide 
an additional method of detection for neutrinos. However, as noted in |29t I38j. the angular sensi- 
tivity for cascades is ~ 50°, which does not permit tracking of the Sun with the desired accuracy. 
Consequently, we do not consider cascade events as a potential signal. 

5 Conclusions 

We investigated the prospects for indirect detection of fermion WIMPless DM at IceCube and 
DeepCore. We considered a hidden sector Majorana DM particle of mass mx that couples through 
Yukawa couplings in the superpotential to a connector of mass my and visible sector particles. 
These models are especially interesting in the context of IceCube/DeepCore because they exhibit 
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only SD nuclear scattering, for which IceCube/DeepCore is expected to soon provide the greatest 
experimental sensitivity for mx ~ C(few 100 GeV). 

We focused our attention on DM annihilation to taus, staus and sneutrinos. In order to be 
captured by the Sun the DM particle needs to couple also to up and down quarks and first generation 
squarks through the superfield. Annihilation to light fermions does not produce an observable signal 
because it is chirality suppressed. As for annihilation into squarks, in the GMSB-inspired WIMPless 
scenario, squarks are generally very heavy and the resulting neutrino spectra are strongly dependent 
on the features of the sparticle spectrum. Therefore, we assumed that the dominant channels for 
annihilation are exclusively leptonic. We assumed that if the DM candidate annihilates to SM 
particles, these are predominantly taus. If it annihilates to supersymmetric particles, these belong 
to lepton superfields and are the NNLSP, with the lightest neutralino being the NLSP. In the cases 
of stau and sneutrino annihilation channels we focused our attention on the distinct situations 
where either the DM candidate couples predominantly to the third family (staus), or it couples 
with equal strength to degenerate sneutrinos of all three families. We believe this allows us to 
consider a range of interesting possibilities. 

We propagated the neutrino spectra originated by DM annihilation at the center of the Sun 
through the solar medium. We took into account energy losses due to NC interactions, CC interac- 
tions, and tau regeneration. We also implemented neutrino oscillations in the propagation through 
the Sun. 

In order to calculate the event rate at the IceCube detector, we considered upward and contained 
events, taking into account muon energy losses due to ionization, bremsstrahlung, pair production, 
and photonuclear effects. The technical details of our analyses are provided in the appendices. 

We found that it is not possible to obtain 3<r detection in upward or contained events at IceCube 
for the stau and sneutrino channel for DM masses mx < 150 GeV with 5 yrs of data. Moreover, 
a 300 GeV DM particle would not produce observable events for the sparticle channels, unless the 
Yukawa couplings are large or the connector mass is not more than about 30% heavier than the DM 
mass. Indirect detection of these relatively light DM particles is favored in the tau channel, due to 
the broader neutrino spectra produced in this channel. We quantified the improvement DeepCore 
brings to the detection prospects in all channels, especially for low mx- In particular, even in the 
cases mentioned above, where the steep decrease in flux below the detector energy threshold makes 
signal detection highly unlikely at IceCube, such limitations do not apply to DeepCore. While the 
performance of the IceCube detector varies significantly between different channels and different 
DM masses, thus requiring a case by case analysis, we showed that DeepCore can comfortably 
produce a 5 yr 3a detection for any analyzed channel, without strong dependence on the Yukawa 
couplings or connector masses. We thus find good prospects for probing models of Majorana 
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fermionic WIMPless DM at IceCube, including the DeepCore extension. 
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A Neutrino Propagation 

We describe the procedure involving the propagation of neutrinos from the center of the Sun to the 
Earth, and the detection at IceCube and DeepCore. 

Once the neutrinos are produced at the center of the Sun, they need to be propagated through 
the solar medium, travel to the Earth and be detected at the neutrino telescope. The appropriate 
formalism involves the density matrix for the neutrino spectra in the flavor basis [39j. We call this 
p, and indicate matrices in boldface. 

Neutrinos of energy E v can be propagated from a point r to r + dr inside the Sun by solving 
the Heisenberg equation, 

dp{E u 



- -i[tL{E v ),p{E v )\ + 



dp{E v 



NC dr 



+ 



CC dr 



(17) 



dr dr 

where H is the Hamiltonian for neutrino oscillations in matter, the term indicated by in is the 
injection spectrum at the center of the Sun and the other two terms represent the matter effects 
due to NC- and CC-interactions. The Hamiltonian is 

H = -^Udiag(0,Am^ 1 ,Am| 1 )U t +diag(v / 2G F A r e ,0,0), (18) 

ZEy 

where U is the neutrino mixing matrix as given in [40J, N e is the radius-dependent density of 
electrons inside the Sun, Gf is the Fermi constant, E v is the energy of the incoming neutrino, and 
Am|j and Am^ are the neutrino mass-squared differences, Amf„- = mf — m]- We use the following 
values: Am^ = 8.1 x 1(T 5 eV 2 , Am| x = 2.2 x 1(T 3 eV 2 , 6 12 = 33.2°, <9 23 = 45°, 6 13 = 0. It has 
been shown in |39} W\\ that the oscillation results are not significantly changed if #13 is small but 
nonzero. The density profile of the Sun is shown in Fig. 12 
The injection term is 



dpi 



dr 



(Eu) = S(r)Si~ . (19) 
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Figure 12: Electron density (solid line) and neutron density (dashed line) as a function of the Sun's 
radius. 



It expresses the number of neutrinos and antineutrinos produced in DM-annihilation at the center 
of the Sun as a function of energy. It is diagonal in the flavor basis. 

In considering the NC and CC terms, we introduce the following quantity which depends on 
the deep inelastic scattering (DIS) total cross sections a: 

Tnc{CC)(Ev, E') = N p {r)d\ag[a{np -> ^'(0 + an y)] + N n (r)di&g[a(v in -> v[{l) + any)] , (20) 

where N p {r){= N e (r)) and N n (r) are the proton and neutron densities inside the Sun, plotted in 
Fig. 12, E v is the incoming neutrino energy, E' the outgoing neutrino (charged lepton) energy and 
I labels flavor. Then, the neutral-current term is given by 



dp{E v ) 



dr 



P{E V ) [ EV ^(E u ,El)dE' u + r^(E' u ,E u )p(El)dE' u . 



(21) 

NC 

The CC term is defined in a similar way but it is more complicated for two reasons. First, 
the CC-DIS cross sections are not the same for all flavors of neutrinos. As a matter of fact the v T 
cross sections are suppressed near threshold by the kinematical effects of m T . Second, one needs 
to take into account the effects of tau regeneration that couple the propagation of different flavors 
(different elements of the density matrix) and elements of the neutrino and antineutrino density 
matrices. 

Tau regeneration is an important effect that leads to a reinjection of the neutrinos produced by 
the decay of the taus that are CC created by neutrinos of higher energies. The taus produced by 
energetic neutrinos undergoing CC interactions can decay promptly through various channels, for 
example, t~ — > ^ r + any, r~ — > e~v e v T and r~ —> jjL~u^v T , and similarly for the antiparticles. These 
processes provide additional sources of energetic neutrinos that reenter the flux with lower energies. 
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The probabilities of reinjection are encoded in four functions f UT ^u T (u), fu T ^u T (u), f Vr - H? e „(«), 
fu T ^u e M ( M )> which depend on the branching ratios of the above channels and on u = E„/E™, 
where E™ is the energy of the tau neutrino undergoing CC scattering, and E™ 1 is the energy of 
the lower energy neutrinos produced by tau decay |39j . The charged-current contribution to the 
Heisenberg equation is therefore, 



dp{E v ) 



dr 



dp(E u 



{Tccp} 



+ 



CC 



00 dE™ 



E' 

— 7 / 



-i- n, l ,p TT (E' l '; )V ; r , .( E'" )f„ T ^ Vailt ( 



dr 



{Tcc,p} 



+ 



CC 



00 dE'J 1 



Win 
E u &v 



U T p TT (E t u n )Tl; C (E l J l )fp T ^p T 



E» 
Ejj 



E v 

+ n eiM p rT (£ , ™)r^ c (^ n )/ I/r ^ PeM ( 



(22) 



(23) 



where U e = diag(l,0,0) are projectors, and similar expressions apply to the other flavors. As is 
clear from the last term on the right-hand side of Eqs. (|22p and (|23p . tau-regeneration effects couple 
the two sets of equations. The regeneration probability functions are given by |41j . 
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(24) 

(25) 

(26) 
(27) 



where z = E T /E™ and N normalizes the equations so that their integral is either 1 {y T 



and v T —7- v T \ or 0.18 {y r 



v e ^ and v T 



'e,n) 



The functions goi an d gu are, respectively, 



the unpolarized and polarized energy spectra of the r neutrinos originating from the taus in the 
fragmentation frame, for each final state. The functions go and g\ are the unpolarized and polarized 
energy spectra of the v e ^ from decay of the r. The explicit forms of goi and gu are given in Table I 
of |43j . and those of go and g\ are given in Ref. [33] • The branching fractions Brj refer to the 

• TT~V T 



six possible final states: r 



e v e v T (Bri =0.18), r — > p v^v T (Br2 =0.18), r 



(Br 3 =0.12), t~ -> oi»/ T (Br 4 =0.13), t" -> pi/ T (Br 5 =0.26), r~ -> i/ r + any (Br 6 =0.13). 

Once the neutrinos reach the surface of the Sun, the propagation to the Earth is obtained by the 
following averaging procedure: we rotate the density matrix to the mass basis; drop the off-diagonal 
terms and rotate it back to the flavor basis. With our choice of neutrino parameters, the averaging 
should wash out any observable modulation. 
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B Muon Rates 



Upward Events at IceCube 

When a muon generated via CC interactions travels through the rock and ice beneath the detector 
it loses energy due to ionization, bremsstrahlung, pair production, and photonuclear effects |45j. 
The average energy loss of the muons that travel a distance dz in a medium of density p me d is given 
by: 

dE\ 

-^)=-(a + P(E)E)p med (z), (28) 

where a = 3.0 x 10 -3 GeV cm /g is related to ionization, while (3(E) takes into account bremsstrahlung, 
pair production, and photonuclear effects. We take j3 = 3.0 x 10~ 6 cm 2 /g and p m ed = Pice = 
0.92 g/cm 3 . The results from Muon Monte Carlo [IB] are reproduced by choosing these values of 
a and (3 [47J. Equation (|28|) can then be easily solved to obtain the final energy E^f, given initial 
energy E^: 

E»f = --p+e-^ z {E, l + fj . (29) 
The average range covered by the muon between energies E^, E„f is then, 

^•*w=sb ta (^)- (30) 

Thus, the muon flux at the detector is obtained by a convolution of the following: the probability 
of the incoming neutrino to CC scatter with a nucleus in ice; the average range over which energy 
losses force the muon energy below the detector threshold; the muon probability of surviving its 
own decay length. This last effect can be parametrized by the survival probability, 

m/j, 

p r -p pl lE^a + PE^* 

^sur y-^/ii > J - J fj.f ) 



CTa Pice 

(31) 



which is a solution to the differential equation, 

dPsur 



(32) 



dE^ f Efj,fCTp ice (a + f3E^f)/m^ ' 
where r is the muon lifetime and its mass. Folding these effects together gives the spectrum of 



muon events, 



^ " / dze^PsuriE^E^^)^^ 
Jo 

M{^ {Eu)Pp + ^E- {Eu)Pn ) 
+ (a corresponding contribution for v)\ , (33) 



dE, 



m x 

dE u 
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where E^E^f , z) is obtained by inverting Eq. ([29]) . da cc /dE u (E u ) are the differential CC-cross 
sections to protons {up) and neutrons (vn), p p ~ 5/9 Na cm -3 and p n ~ 4/9 cm -3 are the 
number densities of nucleons in ice expressed in terms of Avogadro's number N^, and d<& u /dE u is 
the neutrino spectrum at Earth, depicted in our specific cases in the third column of Figs. 3, 4 and 5. 

The event rate for upward events is obtained by convolving Eq. (|33|) with the muon effective 
area of the detector, A e ff(E )1 f)R(cos9), which is constituted by a zenith angle-independent part, 
shown in Fig. 13a, and a factor i?(cos 9) = 0.92 — 0.45 cos 9 that accounts for the rock bed beneath 
the ice |48] . We take the average of the effective area over the time of the year that the Sun spends 
below the horizon, namely between the March and September equinoxes. We define the zenith 
angle 9 Z at the South Pole to be the angle centered at the detector with 9 Z = 0° indicating the 
vertical direction in the sky. 9 Z can be parametrized in terms of the time of the year f y (where 
f y = 0, 1/2 correspond to the March and September equinoxes, respectively), and the tilt of the 
Earth axis with respect to the perpendicular to the ecliptic plane, 9 t = 23° 26': 



Ufy 



7T 



+ 9 t sin(27r f y 



(34) 



The event rate reads, 



JV events ~ <, L eq 



(A eff {E uf )R(cos9 z )}dE af , 



(35) 



'0 dE uf 

where £ and T eq are the quantities introduced in Sec. 3.1. (A e ff(E u f)R(cos9 z )} is the average over 
the portion of the solid angle (9 Z ,4>) that corresponds to the time the Sun spends below the horizon. 
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Figure 13: a) The IceCube muon effective area with zenith angle dependence factored out. b) Muon 
effective volume for DeepCore. 
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Contained Events at IceCube 

For neutrinos that interact within the detector volume, the resulting muons do not lose an appre- 
ciable fraction of their energies. Following notation similar to the previous subsection (except that 
we identify = E^i since the muons do not propagate) the muon flux for contained events is 
given by 

rl& r m X T^CE) / dn CC rlrrCC \ 

(36) 



d ^- L r d E, 



dE n Jeu. 



where L ~ 1 km is the size of the IceCube detector. The event rate reads 



Ng ents = HT eq I J(l km 2 )^ , (37) 



where Ethr = 100 GeV is the energy threshold of the IceCube detector. As mentioned above, we 
only consider events observed between the March and September equinoxes. 

Events at DeepCore 

The great advantage of DeepCore with respect to IceCube is that the outer instrumented volume 
of IceCube will serve as a veto to atmospheric muon events up to one part in 10 6 [29|, so that 
data can be collected throughout the year, i.e., even when the Sun is above the horizon. The 
rate for contained events at DeepCore can be calculated by convolving Eq. (|36p with the muon 
effective volume V e ff(E^i). In the most optimistic estimates the effective volume is constant for 
muon energies above ~ 300 GeV, and drops significantly at lower energies |31j . As explained in 
Sec. 3.3, for DeepCore we consider the interval > E m i n = 35 GeV. We find that in this interval 
the effective volume in km 3 can be parametrized by 

V eff (E^) = (0.0056 log E„i + 0.0146)6(275 - E^) + 0.02839(£^ - 275) , (38) 

where is the Heaviside step function and E^ is in GeV. The effective volume is plotted in Fig. 13b. 
After convolution one gets 

Kents = / jj^VeffiE^dE, . (39) 

C Atmospheric Background 

The angle-dependent flux of atmospheric neutrinos d^ tm /(dE u d cos 9) is given in [33 . These 
neutrinos interact with the medium surrounding the detector and produce muons that constitute 
the background. As in the case of neutrinos from annihilation, the atmospheric background can be 
divided in upward and contained events. 
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By following the notation introduced in the previous subsections we can write the angular 
dependence of the upward background flux for IceCube as 

I 1 



dE^fd cos 9 z (f, 



/ dze^^PsuriEiiE^.z)^^ 
Jo 

Emcix 



d$l tm ( d<jCC drfc 
dEyd cos (j z \Jy) \ dEi dEi 



dE v 

+ -> v)\ , (40) 

where Emux ~ 10 4 GeV is the maximum energy above which the rapidly falling atmospheric flux 
becomes negligible. The angle-dependent flux, Eq. (|40p . is integrated over a cone of solid angle 
dJ7 = ir(d9) 2 , where the opening angle (half the apex angle) d9 ~ 1°, consistent with the IceCube 
angular sensitivity [29(, [M] . The cone tracks the Sun according to Eq. (|34p , for the fraction of the 
year over which the Sun is below the horizon. The number of upward background events for IceCube 
can be calculated by folding in the effective area A e ff(E /1 f)R(cos9 z ). We find N B q ~ 6.1 yr -1 
where, again, data is taken only between the March and September Equinoxes. 

Contained events are those obtained by muon neutrinos undergoing charged-current interactions 
within the detector volume, thus the angular dependence of the flux is given by 



- L I dE u ( p p — + p n — + 



atm 
v 



dE ^d cos 9 z {f y ) J Efi V dEfj, ' dE^ J dE v d cos 9 z (f y ) 

+ {v^v). (41) 

The number of contained background events in a cone of 1° opening for time of exposure limited to 
the period between March and September is obtained, as in Eq. (|37p . by convolving with a constant 
1 km 2 area. We get N^q 1 ~ 15.6 yr -1 . 

The number of background events for DeepCore over the whole year, obtained by convolving 
with V eff (E^)/L, is Ngg ~ 2.5 yr- 1 . 

In Fig. 14 we show the dependence of N^q, N^q 1 an d ^bg on ^ ne °P enrn g angle d9. The thin 
blue horizontal lines shown in Figs. 6-11 represent the 5 yr 3<r-discovery reach for d6 = 1°. While 
the angular resolution is expected to be best for directions close to the horizon (and the Sun does 
not stray more than 23° 26' from the horizon at the South Pole), in case the angular resolution of 
the IceCube and DeepCore detectors do not meet expectations, the background event rates can be 
obtained from Fig. 14. Then, for a 3<r detection, iV„ = 3\/Nbg can be trivially recalculated. 
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Figure 14: Number of background events per year as a function of the opening angle d6. 
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